Local Authority Boundary with Air Quality Management Areas
leaflet() %>%
addTiles() %>% # Add default OpenStreetMap map tiles
addPolygons(
data=AQMA[zones[i,],], weight = 1) %>%
addPolygons(
data=zones[i,], fillOpacity = 0, weight = 2)
© Ordnance Survey & Crown copyright and database right 2017
NOx
if(!is.null(nox)){
plot(nox$date, nox$nox)
} else {
print("No NOx measures")
}

# df <- get1Hdata('MY1', years=2015)
#
# # Aggregate to daily means and plot
# library('zoo')
# par(mai = c(0.5, 1, 0, 0))
# my1 <- zoo(x = df$Ozone, order.by = as.POSIXlt(df$datetime))
# plot(aggregate(my1, as.Date(as.POSIXlt(df$datetime)), mean),
# main = '', xlab = '', ylab = expression(paste('Ozone concentration [',
# mu, 'g/', m^3, ']')))